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Abstract 

We solve numerically a field equation of axions coupled with gravity and 
show solutions representing oscillating axionic boson stars with small mass 
~ 1O~ 12 M0 and large radius ~ 10 8 cm. We present explicitly a relation 
between the mass and the radius of the boson stars with such a small mass. 
These axionic boson stars are shown to possess oscillating electric currents 
in an external magnetic field and to radiate electromagnetic fields with a 
frequency given by mass of the axions. We estimate roughly the luminosity 
of the radiation in a strong magnetic field 10 12 G of neutron stars, it is ~ 10 34 
erg/s with the axion mass 10 -5 eV. 



The axion is a Goldstone boson associated with Peccei-Quinn symmetry |T| and is one of 
most plausible candidates of the dark matter in the Universe. The axions are produced 
|||| in early Universe mainly by the decay of the axion strings or coherent oscillation of the 
axion field. These axions form incoherent axion gas in the present Universe. It has been 
recently argued [|]] that not only the incoherent axion gas but also coherent axionic boson 
stars are produced by the coherent oscillation of the axion field. Namely, inhomogeneity 
of the coherent oscillation on the scale beyond the horizon generates axion clumps as well 
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as axion domain walls. The domain walls decay soon after their productions and leave a 
primordial magnetic field || as well as incoherent axions. On the other hand, the axion 
clumps contract gravitationally in later stage of the Universe to the axionic boson stars. We 
call them simply as axion stars. 



of the axion together with Einstein equations. Our concern is confirming the existence of 
the axion stars with quite small masses and finding the explicit relation between the radius 
( or the typical value of the axion field ) and the mass of the star. This is because the mass 
of the axion star produced by the mechanism mentioned above has been estimated to be 
order of 1O~ 12 M by Kolb and Tkachev Furthermore, using these solutions, we show 
that the axion stars in an external magnetic field can emit electromagnetic radiations with 
a frequency given by the mass of the axion. The luminosity of the radiations is roughly 10 34 
erg/s in the magnetic field ~ 10 12 G of neutron stars. 

Now we show solutions of the axion stars by analyzing the field equations, 



where we have assumed gravity being small, i.e. h t>r C 1 so that the metric is such that 
ds 2 = (1 + h t )dt 2 — (1 + h r )dr 2 — r 2 (d9 2 + sin 2 9d(p 2 ); r,9, and 4> denote the polar coordinates. 
The first equation is the equation of the axion field a coupled only with gravity. The second 
and the third equations are Einstein equations. A dot ( dash ) indicates a derivative in time 
t ( r ). G ( m ) is the graviational constant ( the mass of the axion ). We have neglected the 
potential term of the axion because an amplitude of the field a is assumed to be sufficiently 
small for the nonlinearity of a not to arise. This implies that the mass of the axion star is 
small enough. Actually the masses we are concerned with are such as ~ 10~ 12 M Q . We need 
to impose a boundary condition such as h r (r — 0) = for the regularity of the space-time. 
Changing the scales such that r = mt, x = mr and b = a/m, we rewrite the equations 



In this letter we present numerical solutions of the axion star by solving a field equation 
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as follows, 



b = Vb + b" + (- + V')b' -b , (4) 

x 



V = K - K = c{ l0 d ^ 2 ( b ' 2 + b2 + b2 ) _ xb 2 } (5) 



with e = AnGm 2 , where we have expressed V in terms of the field b, solving eq(§) and eq(^); 
here a dot ( dash ) denote a derivative in r ( x ). We understand that if the gravitational 
effect is neglected ( e = ), the equation of b is reduced to a Klein-Gordon equation. Thus 
the frequency uj of the field b is modified from the value of the free field, 1 ( m in the physical 
unit ), only by a small quantity proportional to the gravitational effect, e; uj = 1 — o(e). We 
look for such a solution |7| that 

b = A B(x) sin ujt + o(e) sin 3ujt , (6) 

where B(x) represents coherent axions bounded gravitationally with its spatial extension 
representing the radius of the axion star; B(x) is normalized such as B(x — 0) = 1. Later 
we find that A is a free parameter determining a mass or a radius of the axion star. The 
second term is a small correction of the order of e. Inserting the formula eq(^) into eq([|) 
and eq(||) and taking account of the gravitational effects only with the order of e, we find 
that 

k>B = B» + (- + «TAB' + ^ ( ! } ^ B (7) 

with 

T^^4^ and U= dy ^ ZB 2 ™ - yB 2 ) , (8) 
x 2 Jo y 2 

where k 2 ( = 1 — uj 2 ) is a binding energy of axions. We have imposed a boundary condition 
for the consistency such that V(x = 0) = h t (x = 0)/2 = eA^v uj sin 2ujt ; this is the definition 
of constant v in the above formula. 

We can see that the parameter eA\ can take an arbitrary value and that it represents the 
gravitational effect of this system. Namely the mass of the axionic boson star is determined 



by choosing a value of the parameter. Note that the normalization of B has been fixed in 
eq([?P although the equation is a linear in B. 

We may take the value of v without loss of generality such that v = —U(x = oo). Then 
the inverse k~ x of the binding energy is turned out to represent a radius of the axion star; B 
decays exponentially such as exp(— kx) for x — > oo. It turns out from eq(J7|) that the choice 
of small values of eA\ lead to solutions representing the axion stars with small mass and 
large radius k^ 1 . 

Before solving eq(0) numerically, it is interesting to rewrite the equation as following. 
That is, we rewrite the equation by taking only dominant terms of a "potential", Vb, in 
eq®, 

^(<T + f, B ' + <^, (9) 

in the limit of the large length scale; setting x = Xy, we take a dominant term as A — > oo. 
This corresponds to thinking the axion stars with their spatial extension being large. Note 
that our concerns are such axion stars with small mass and with large radius. 

Then, since the dominant term in the limit is the last term in eq(|S]), U + v ~ J£° dxxB 2 , 
we obtain the following equation, 

- -„ 2B' BJ°°dyyB 2 . . 

B = B" + + — ± — — — , (10) 

z 2 

where we have scaled the variables such that B 2 = k 2 B 2 jeA\ and x = k~ x z; a dash denotes 
a derivative in z. This equation is much simpler than eq(|7]), where we need to find each 
eigenvalue of k for each value of eA\ given, in order to obtain solutions of the axion stars 
with various masses. On the other hand, we need only to find an appropriate value of 
B(z = 0) = k 2 /eAQ in order to obtain such solutions. A relevant solution we need to 
find is the solution without any nodes. Obviously, the solution is characterized by one free 
parameter, k 2 or eA^, which is related to the mass of the axion star. Namely the choice of 
a value of the mass determines uniquely the properties of the axion star, e.g. radius of the 
star, distribution of axion field a, e.t.o. 



Although the equation (|T0|) governs the axion star only with the large radius, we have 
confirmed by solving numerically original equation (0) that the stars of our concern, whose 
masses are at most 1O _12 M , can be controlled by the equation (|T0|). Thus we may determine 
the explicit relation between the mass and the radius of the axion star, by finding the above 
solution of eq(]T0|) , 
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M = 6A^- , (11) 
m 2 R 

where m v \ is Planck mass. Numerically we can see that for example, R = 1.6 x 10 5 m^" 2 cm 
for M = 1(T 9 M , R = l.Qx 10 8 m5 2 cm for M = 1(T 12 M , e.t.c. with m 5 = m/l(T 5 eV. 
We have also numerically confirmed that the axion field a of the solution can be practically 
parameterized such that a = f a exp(— r / R) where / is the decay constant of the axion 
whose value is constrained from cosmological and astrophysical considerations 0; 10 10 GeV 
< / < 10 12 GeV. This implies a constraint on the mass of the axion; 10 _5 eV < m < 10~ 3 eV. 
Amplitude a ( = mA / f ) is found to be expressed in terms of the radius R, 

s (10 8 cm) 2 10" 5 eV 

ao = 1.73x 1(T 8 1 ; . 12 

R z m 

In this way we have found solutions representing the axion stars with their mass M being 
typically 1O~ 12 M and with the radius R being 10 8 cm. The binding energy is turned out to 
be extremely small so that the frequency u of the field a is given by m. 

Next, we wish to show that the axion stars can emit electromagnetic radiations with 
energy m in an external magnetic field. The reason is that the coherent axion field gain an 
oscillating current in the magnetic field owing to the following interaction, 

caaE ■ B 

with a = 1/137, where E and B denote electric and magnetic field, respectively. The value 
of c depends on the axion models |S|||; typically it is the order of one. 

It is easy to see from this interaction that the coherent axion field may have electric charge 



density —cada ■ B / fn in the magnetic field B [ 1Q| . We assume that the field B = (0, 0, B) 



is spatially uniform and that the geometry of the axion field a representing the boson stars 
is spherical. Then we understand easily that this axion field has a charge distribution such 
that it has negative charges on a hemisphere ( z > ) and positive charges on the other 
side of the sphere ( z < ). Net charge is zero. Thus, the star possesses the electric field, 
E, parallel to the magnetic field associated with the charge distribution. As we have shown 
before, the field a representing the axion star oscillates with frequency being approximately 
given by m. Therefore, the charge distribution oscillates with the frequency and so an 
electric current associated with the charge oscillates similarly. This fact induces emission of 
the electromagnetic fields. 

The radiation from the axion stars can be understood in another way. Since a mix- 
ing angle of an axion and a photon generated by the interaction eq(|13D is the order 



of caB / fmir in the magnetic field B, equations of motion of gauge potentials Aj can be 
derived such as d 2 Ai = d 2 caBia/ fmir ~ caBimaj fir. Hence, the oscillation of this current, 
caBima/ fn, produces the radiation as dictated above. 

It is easy to evaluate the luminosity L of the radiation by assuming the configuration of 
the field a such that a = faoexp(—r/R) whose form we have confirmed numerically. Noting 
the wave length ( = m _1 ) of the radiation is much smaller than the radius R of the axion 
star, we obtain 

L ~ 3m*R 2 ' ( } 

where we have assumed that whole of the axion star is involved in the external magnetic 
field. On the other hand, when the length scale of the magnetized objects like neutron stars 
is smaller than the axion stars, R in the above formula eq(|14D should be understood to 
denote the radius of the magnetized objects. Especially what we are concerned with is the 
collision of a neutron star with the axion star whose radius is 10 8 cm ~ 10 10 cm. This radius 
is 10 2 ~ 10 4 times larger than the radius of the neutron star. Thus R in eq(0) denotes the 
radius of the neutron star in such an example. 

We understand that since the luminosity L is proportional to B 2 , it is larger as the axion 
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star is exposed to a stronger magnetic field. Furthermore, as L depends on R~ 2 , it is larger 
as the radius of the axion star or the magnetized object is smaller. Hence we understand 
that the neutron star among astrophysical objects yields large luminosity when it collides 
with the axion star. We estimate the luminosity of the case, 

, 9 , B 2 M 4 (10 6 cm) 2 , . 

L = 2.7xl0Verg/s (1|)12G)2(10 _ 12Ms) / ff , (15) 

where we have used the above formulae egflTl]) and eq(|l2D for expressing ao in terms of 
the mass M of the axion star. We note that this luminosity is that of the monochromatic 
radiation with the frequency, m/2ir = 2.4 x 10 9 Hz (m/10~ 5 eV). This is much weaker than 
a corresponding luminosity of synchrotron radiations with the frequency from the neutron 
stars. This emission continues until the axion star passes through the neutron stars. It takes 
10 8 cm x(lO~ 12 M /M)/(3 x 10 7 cm/s) ~ 10 sec (1O~ 12 M /M), assuming that the velocity 
of the axion star is 3 x 10 7 cm/s, which is the typical velocity of matters composing the halo 
in our galaxy. Thus it is impossible to observe the radiation. 

Until now we do not take into account the influence of the medium of the astrophysical 
object such as the neutron star. When the medium is electric conducting, we need to consider 
its effects on the radiations or the electric field associated with the axion star. Here we would 
like to mention that the electric field possessed by the axion star in an external magnetic field 
induces an electric current in an astrophysical medium. This current also oscillates and emits 
a radiation. Since the strength of the current is proportional to an electric conductivity, a, 
of the medium, the medium with large conductivity such as the crust of the neutron star 
leads to strong radiations. In the case, the luminosity is approximately given by 

L a = (-) 2 L (16) 
m 

where the conductivity should be taken as an conductivity averaged over the volume of the 
medium. For instance, in the case of the neutron star the conductivity of its crust has been 
estimated [JEJ theoretically such as a =O(10 24 /s). When we adopt this value for a, the 
luminosity L a is given by ~ 10 34 erg/s. Although this is a rough estimation, we can expect 



that the neutron star emit sufficiently strong radiations to be detectable when it collides 
with the axion star. The detail analysis is in progress. 

Finally we point out that the monochromatic radiation under the discussion should be 
detected for confirmation of the phenomena associated with the axion just after the discov- 
ery of rebrightness of dark neutron stars or white dwarfs caused by the axion star. The 
rebrightness of these stars has been discussed in the previous paper |13|| to be caused by 
energy dissipation of the axion star when it collide with them. These two phenomena, 
monochromatic radiations and thermal radiations owing to the energy dissipation, are dis- 
tinctive in the collision of the axion star with these strongly magnetized stars. 

In conclusion, we have shown the solutions of the oscillating axionic boson stars with 
small masses. We have made explicit the relation between the mass and the radius of these 
axion stars. They have also been shown to gain two types of oscillating electric currents in 
the magnetized conducting medium; currents made of axions themselves and currents made 
of electrons in the medium. Both of them emit monochromatic radiations with the frequency 
given by the mass of the axion. Having estimated the intensity of the radiations, we found 
that the radiations by the currents of electrons are strong enough to be detectable. But, 
since the duration of the emission is much short, it seems difficult to observe the radiations. 
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